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One can describe cosmological relic neutrinos by adding Lagrange multipliers to the Standard 
Model Lagrangian for them. The two possible Lagrange multipliers are a chemical potential, which 
fixes the mean neutrino/anti-neutrino asymmetry, and a Majorana mass, which fixes the mean spin- 
entropy. Because these neutrinos originated from a thermal bath, their entropy should be maximal, 
implying that each state in the background is a symmetric superposition of a neutrino and anti- 
neutrino. Therefore the Standard Model must be augmented by a flavor-diagonal Majorana neutrino 
mass matrix. This impacts the propagator via tadpole diagrams due to self-interactions. In the low- 
energy limit, neutrino self-interactions are entirely off-diagonal because same- flavor four-fermion 
operators vanish by Pauli exclusion. These interactions must be diagonalized when propagating 
through a bath of neutrinos, using the J7(3) global flavor symmetry. U{?>) gets broken broken 
down to SO(3) by Majorana masses, and down to A4, if the three masses are different. Thus our 
universe today contains tri-bimaximal mixing and Majorana neutrinos. Neutrino mixing is due to 
the mismatch between the flavor-diagonal Majorana mass matrix arising at finite density and the 
self-interaction diagonal finite density propagator. The mass hierarchy is inverted and Majorana 
phases are absent. Lepton number is conserved and the neutrino-less double beta decay experiment 
absorbs a pair of neutrinos from the relic background and will prove their Majorana nature. 



INTRODUCTION 

Neutrinos have been observed to change flavor. Origi- 
nally known as the "solar neutrino problem" due to the 
apparent deficit of solar neutrinos, today this is well- 
understood as the mixing of electron-type neutrinos emit- 
ted in nuclear reactions into the /x and r types. The 
most economical way to describe this mixing is via an 
off-diagonal mass matrix. One can implement this with 
either Majorana or Dirac neutrino masses. We currently 
have no experimental evidence that neutrino masses are 
either Dirac or Majorana, though there are several run- 
ning and planned experiments to test this using neutrino- 
less double beta decay, which can only occur with a Majo- 
rana mass.flj The current experimental situation is well- 
explained by a three-neutrino mixing hypothesis. 

The standard cosmological model predicts the exis- 
tence of a sea of relic neutrinos if the temperature was 
ever above a few MeV. The leading order dynamics of this 
background that is relevant for neutrino experiments is 
its Pauli blocking. This can be included by adding La- 
grange multipliers to the Lagrangian that fix the number 
density and entropy. 

Whenever particles propagate in a medium, an index 
of refraction and rescaling of the kinetic term are in- 
duced, defined by the tadpole diagram corrections to the 
self-energy. This is a non-perturbative effect that must 
be re-summed using the Dyson series into the two-point 
function. 



MEAN FIELD NEUTRINOS 

Our universe contains neutrinos, therefore when doing 
"vacuum" physics we should fix the parameters of the 



background relics with which we arc not scattering. This 
is the approach taken whenever one adds a chemical po- 
tential directly to the Lagrangian. One is averaging over 
the degrees of freedom represented by the chemical po- 
tential, treating them instead as a "mean" field. One can 
also envision this as constraining the first moment of the 
thermal distribution function. 

A finite density system in an arbitrary configuration 
can be represented by a density matrix 



(1) 



The indices i,j represent the distinct quantum numbers 
possible for a single momentum mode p, such as helicity 
and flavor. The average value of this density matrix is 



{P^. 



(2) 



In the mean field approach one assumes that the back- 
ground state pij is independent of time. By then adding 
Lagrange multipliers for the components of pij (j>) (writ- 
ten in terms of the field operator) , one ensures that the 
result Eql2] is reflected in the dynamics. This assumes 
that Eql2]will be unchanged by the dynamics we are in- 
terested in computing. 

Let us consider a single massless Weyl fermion. Its 
2x2 spin-density matrix operator p can be related to 
the field operators using 



{a+, a-); 
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where a± are the annihilation operators for the two helici- 
ties, and the momentum dependence has been suppressed 
for notational simplicity. The field operators for a Weyl 
fermion are 



-ipx 
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with the commuting state vectors in the helicity basis for 
neutrinos and anti-neutrinos respectively 

In terms of the density operator Eql3l wc can find the 
relationship between the density matrix and Lagrangian 
operators by evaluating Tr[/5cr-']: 

Tr[paO] = I 



Tr[/5ai] 

Tr[pa2] 
Tr[pa3] 



(6) 



The first of these specifies the normalization of the den- 
sity matrix. The second and third combinations arc both 
Majorana masses, and are related by a Icpton number 
phase transformation. Finally the fourth is the number 
operator. 

In the mean field approximation described by the 
momentum-averaged density matrix, Eql^l one can now 
add ^jTi[paj] to the Lagrangian with constant coeffi- 
cients. This is a Lagrange multiplier, and when finding 
the equations of motion for the fermion Xi will ensure 
that these average values specified by the density ma- 
trix are preserved. Note that does not transform as a 
Lorentz vector. 

Let us examine further the nature of these operators. 
With two neutrinos in a volume V, the averaged density 
matrix is 



(Ph 



1-Ri J V 
The appropriate traces analogous to Eql6] are 



(7) 



Tr[(p)2ai] = 2i?2^1-i?2cos(^2)_ (g) 



Tr[(p)2a2] 



2i?2A/l-i?isin(^2)-^ 
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(9) 
(10) 



where we've taken a unit volume V. For a pure 
state of two neutrinos, i?2 = 0, the Majorana mass 
term Tt[(p)2(J^] = 0, and the chemical potential term 
Tr[(p)20''^] = —2/V is negative. For a pure state of 
anti-neutrinos, i?2 = 1 and the chemical potential term 
changes sign, while the Majorana mass term remains 
zero. For a mixed state with an equal probability for 
each to be a neutrino or anti-neutrino, i?2 = 1/V2, the 
Majorana mass terms are maximal, and the chemical po- 
tential term Tr[(p)20''^] is zero. 

The chemical potential /is is expected to be related to 
the baryon asymmetry of the universe, and very small 



6 X 10"^". 



Therefore we will neglect it and assume /is = so that 
the number of neutrinos and anti-neutrinos are the same. 

The state with both /ti and /i2 zero is a superposition 
of a pure state of neutrinos with a pure state of anti- 
neutrinos. The averaged density matrix is proportional 
to the identity. This configuration has zero spin-entropy 



Si 



-kTr[{p)ln{p)]. 



(11) 



Since the early-universe plasma from which these neutri- 
nos originated should be extremely close to a blackbody 
spectrum, as is observed for the Cosmic Microwave Back- 
ground, the entropy for this system should instead be 
maximal. For the spin-entropy this corresponds to max- 
imizing the magnitude of the off-diagonal components, 
e.g. R ~ I/V2, and is a state that is in an equal super- 
position of neutrino and anti-neutrino. 

If we examine the thermal Fermi-Dirac distribution 
f{E) = (e^/'^"^ -I- for a system with zero chemi- 
cal potential, we see that f{E) < 1/2 for all energies. 
In other words, each momentum mode there is either a 
neutrino or an anti-neutrino, but not both (which would 
correspond to f{E) = 1 for that energy). If both were 
present for a single momentum mode, Fermi statistics 
indicates that the two density matrices are orthogonal 
Pi ■ P2 = 0, which implies that their sum is the identity. 
This configuration has zero spin-entropy also. 

Each pair in the symmetric bath therefore contributes 
S± = A: 2 In 2 to the spin-entropy. Expanding in Eq llll 
(p) = 1 + e, where e is off diagonal, we can see that only 
the magnitude of the off diagonal components of (p) con- 
tribute to the spin-entropy. If the phase were random, 
as in the Random Phase Approximation (RPA) , these off 
diagonal elements would cancel in the sum, Eq[2] leading 
again to zero spin-entropy. This is clearly incorrect. In- 
stead the condition of maximal spin-entropy implies that 
all states have a fixed, common phase. 

The two masses pi and p2 are related by the lepton 
number phase symmetry x — > e*''x- The Standard Model 
conserves lepton number for each species individually, so 
we can use this freedom to rotate away one of pi or p2. 
In other words, the Majorana phases a = arctan(/ii//i2) 
if present, can be set to zero and are unphysical. Only 
Ml + is physical. 

Thus given the expectation that the relic neutrino 
background has small chemical potential and maximal 
entropy, we should add to the Standard Model a Ma- 
jorana mass (which we can choose to be pi or P2) for 
each fiavor, as well as a small chemical potential. We do 
not add flavor off-diagonal terms because the creation 
and annihilation operators in Eql6] are related by the 
Lorentz representation of the Weyl field. In other words, 
a Lorentz transformation docs not perform flavor rota- 
tions. 

Thus, the Majorana mass should not be considered a 
"fundamental" Lagrangian parameter. Instead it is an 
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environmental parameter; a Lagrange multiplier analo- 
gous to the chemical potential for a system with spin- 
entropy. Because the spin-entropy is maximal and sim- 
ply related to the number density, the Majorana mass is 
also in effect fixing the total number of neutrinos. Quan- 
titatively then, it is simply the Fermi momentum 



(12) 



where n is the number density. 

SELF-INTERACTIONS 



Neutrino self-interactions occur at tree level by Z bo- 
son exchange. The energies of experiments which have 
demonstrated neutrino mixing arc far below the energy 
where the process vv —f vv becomes large because the 
Z boson is on shell. Therefore we may integrate the Z 
out using its leading order equation of motion at zero 
momentum, which generates four-point neutrino contact 
operators. 

For neutrino energies E <C M'^/T^, only these interac- 
tion terms survive 

Jr{Ji)x = (xk^Xe + xXXp + xt^^Xr)' (13) 

in terms of the currents J/^ = e*"'''x]^^Xfc- The same- 
flavor four-fermion interactions vanished in this equation 
because 



xV^XX^c^AX 







(14) 



which is another way of saying that same-flavor interac- 
tions only occur in the p-wave due to Pauli blocking, and 
vanish in the low energy limit below the Z mass. 

The interactions of Eq[T3] clearly do not correspond to 
the flavor-diagonal kinetic terms. Therefore we should do 
a basis rotation to make them flavor diagonal before com- 
puting the Green's functions, using the Standard Model's 
vacuum global symmetry C/(3) on the three families of 
leptons. Adding Majorana masses that are the same for 
all families breaks all the real generators of this symme- 
try. The remaining pure imaginary generators form a 
representation of 5*0(3). 

This residual 5*0(3) of the Standard Model with three 
equal Majorana masses and its Noether current are 



Xi 



e-^x,; 



Jt = eyfex]^\fe, (15) 



where the Latin indices run over the three families 

and the three generators are 
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- Xl-^\e 




= xl^^Xr 


- X^^^Xm 


■-'Te 


= xl'^^Xe 


- Xk\r. 



with a = {efi, fjiTjTe}, and we can see that the currents 
in Eq[T5]are off-diagonal relative to the flavor basis. 



(16) 



Using our vacuum flavor 5*0(3) freedom, we can find a 
basis in which the interacting currents correspond to the 
kinetic term fields (which remain diagonal under a 50(3) 
flavor rotation). Because the couplings are all the same 
there is nothing to differentiate the flavors in the vacuum, 
so we know that the diagonal interaction will take the 
form of a sum of squares of a current, and that current 
must be the same for all three fields. The rotation that 
accomplishes this is the "tri-bimaximal mixing" matrix, 
and has the form 




U= L J_ J_ 

"75 "73 V2 
The three identical currents after the rotation are 
1 



(17) 



J. 
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(18) 



leading to an effective interaction 



9z_ tA Tint _ 
2 •^int'^X ^ 



2M 
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^xl^^X^xl^AX*. (19) 

i=l 
3 

i^^l^PL^l'^VlxPL^^ = 0. 



where we have written the expression using real Majo- 
rana spinors in the last line. This is the vacuum self- 
interaction basis of neutrinos in the Standard Model. 
However one can see by EqlH] that because different fla- 
vor fermions also normally taken to anti-commute, this 
rotated four-Fermi operator vanishes too! Therefore, at 
0{Gf) neutrinos have no self- interactions. 

Eq |19l has a discrete A4 symmetry, in which the legs 
of the four-Fermi operator are exchanged (this is easier 
to see using real Majorana spinors). The discrete sym- 
metry of the interchange of four objects generates the 
"tri-bimaximal" mixing matrix, EallTI 

The interaction term of EqlT3] and Eq[Tn] arc identical. 
Nothing in the vacuum Lagrangian forces us to choose 
this, or any other basis. It is the basis defined by U 
that the neutrino self-interactions are diagonal in the 
low-energy limit, and the neutrino propagator is most 
transparently defined when propagating through a neu- 
trino medium. One can think of neutrinos as sharing 
one propagator, even though they interact with charged 
currents as three separate fields. 
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Following the arguments of the previous section, we 
must now add a flavor-diagonal Majorana mass matrix 
to the Lagrangian. We chose the flavor basis because 
one should first write the vacuum Lagrangian in a basis 
with diagonal 2-point correlation functions for all fields 
(especially the charged leptons) and then add Lagrange 
multipliers corresponding to the degrees of freedom in 
that vacuum Lagrangian. Writing masses diagonally in 
any other basis is equivalent however. 

This mass matrix does not have the SO (3) symmetry 
of the vacuum Lagrangian, and breaks the 50(3) down 
to the discrete A4 in Eq[T9l If two of the masses were 
equal there would be a remaining S0{2) symmetry, but 
we already know experimentally this is not the case. 



symmetry that we found among the neutral interactions 
at 0{Gf)- Therefore, all three of these generate correc- 
tions to U at 0{Gj^). The largest of these arises from 
the T Icpton, which generates corrections 

The contributions at 0{Gf) found in Refs.[l,0] vanish 
by Eg 1141 in the case of same-flavor interactions and are 
only present when considering a fictitious experiment of 
a neutrino propagating through a bath of neutrinos of a 
different flavor. Radiative corrections to these terms are 
also computed in Ref . . 



NEUTRINO 2-POINT GREEN'S FUNCTION 

The purpose of the analysis in the two preceding sec- 
tions is to enable us to compute the full, in-medium two- 
point correlation function. This means computing the 
one-particle irreducible (IPI) tadpole interaction with 
the background. The effect of these tadpoles is inher- 
ently non-perturbativc, and cannot be neglected due to 
a small coupling with the medium. 

The propagation of a single fermion is governed by the 
equation of motion 

[ft - S(fc)] V = 0. (20) 

The self-energy S(fc) has the general form 

Y.{k)=m-a^ + bi^ + c[^,iL] (21) 

where u" is the frame in which the thermal distribution 
function is deflned.Q Terms proportional to can only 
show up at 0{Gp) in the low energy theory, because the 
Majorana mass is Lorentz invariant. It is an expectation 
value for the energy (E) ^ that defines the rest frame 
of the medium u". Here a contributes to a rescaling of 
the kinetic term, and h and c contribute to the index of 
refraction. 

The really dangerous term that necessitates the diag- 
onalization of Eq[T7| is a. If we had done these com- 
putations in a basis with an off-diagonal mass matrix 
(whatever its origin), and off-diagonal self-interactions 
we would have also generated off-diagonal kinetic terms, 
via the a and c of Eq l21l that must be diagonalizcd in 
order to have a sensible set of correlation functions. 

Using the self-interaction basis in Eq[T7l the neutral 
current interactions are entirely removed at 0[Gf\ The 
finite density contributions to Eq[2T] therefore arise at 
0(G\^ due to charged bosons and p-wave Z bosons. 
We have kept the kinetic term diagonal via this rotation 
at the expense of generating non-diagonal matrices for 
neutrino mass, neutrino chemical potential, and charged 
Icpton mass. All three of these matrices break the Ai^ 



DISCUSSION 

The presence of the relic neutrino background in 
the universe necessitates diagonalizing neutrino self- 
interactions in order to compute the neutrino propaga- 
tor. The dominant observable consequence of this is an 
apparent flavor mixing due to the mismatch between the 
propagator's self-interaction basis and the flavor basis in 
which neutrinos are created today by bosons, and 
their spin-entropy defined. 

The mixing matrix is given uniquely by the require- 
ment that the interactions are diagonal, and is of "tri- 
bimaximal" form, previously guessed but not explained 
by other authors. 0, Si Therefore 

sin^ 6*12 = 1/3 sin^ 6*23 = 1/2 sin^ 6*13 = (23) 

in the Standard Model, in the absence of new physics, 
if there exists a relic neutrino background containing 
all three fiavors. The leading corrections to this are 
0(77i2GF/167r2) ~ 10-7. 

In the mixing matrix sin^ ^13 = also implies the ab- 
sence of CP violation, which can be understood because 
the relic neutrino background with no chemical poten- 
tial is itself C and CP symmetric. A tiny sin^ ^13 will be 
generated by loops of charged leptons, and a further cor- 
rection as well as an extremely tiny CP violating phase 
will be generated by a chemical potential /^3. 

Because the freezes out of thermal equilibrium be- 
fore the and tt^ , the relic density of is smaller than 
that of Vf^ and i'^ today. The tt^ and /i^ are very similar 
in mass, so they freeze out creating similar densities of 
Vf. and v^^. Thus due to cosmic expansion the density of 
Vt has had more time to dilute than that of Vf^ and j^^, 
so the hierarchy is inverted with (5tosoi = — m^^ and 

We unambiguously predict that the Standard Model's 
neutrino mass is of Majorana type. Therefore neutrino- 
less double beta decay experiments will observe the neu- 
trino mass, and can be understood as absorbing a pair of 
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neutrinos from the background. New physics which gen- 
erates mixing is perturbing around tri-bimaximal mixing, 
rather than no mixing. 

ACKNOWLEDGEMENTS 

We thank Jose Bernabeu, Kiwoon Choi, Radovan 
Dermisek, Joachim Kopp, Pasquale Serpcio, Thomas 
Schwetz and Jure Zupan for valuable conversations and 
comments on this manuscript. 



* lbob.mcelrath @cern.ch' 
[1] M. Doi, T. Kotani, and E. Takasugi. Double beta Decay 
and Majorana Neutrino. Prog. Theor. Phys. SuppL, 83:1, 
1985. 

[2] Thomas Schwetz, M. A. Tortola, and Jose W. F. Valle. 
Three-flavour neutrino oscillation update. New J. Phys., 
10:113011, 2008. 



[3] H. Arthur Weldon. Effective fermion masses of order 
gT in high-temperature gauge theories with exact chi- 
ral invariance. Physical Review D, 26(10):2789, November 
1982. Copyright (C) 2009 The American Physical Society; 
Please report any problems to prola@aps.org. 

[4] Dirk Notzold and Georg Raffelt. Neutrino dispersion 
at finite temperature and density. Nuclear Physics B, 
307(4) :924-936, October 1988. 

[5] Juan Carlos D'Olivo, Jose F. Nieves, and Manuel Torres. 
Finite temperature corrections to the effective potential of 
neutrinos in a medium. Phys. Rev., D46:1172-1179, 1992. 

[6] Paul Langacker and Jiang Liu. Standard model contribu- 
tions to the neutrino index of refraction in the early uni- 
verse. hep-ph/9206209, June 1992. Phys.Rev. D46 (1992) 
4140-4160. 

[7] P. F. Harrison, D. H. Perkins, and W. G. Scott. A re- 
determination of the neutrino mass-squared difference in 
tri- maximal mixing with terrestrial matter effects. Physics 
Letters B, 458(1):79 - 92, 1999. 

[8] P. F. Harrison, D. H. Perkins, and W. G. Scott. Tri- 
bimaximal mixing and the neutrino oscillation data. 
Physics Letters B, 530(1-4):167 - 173, 2002. 



